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Abstract 

The asymptotic behaviour at large N of the MS quark anomalous dimensions is derived 
to all orders assuming only MS factorization and standard results for the exponentiation 
of soft logarithms in the quark initiated bare cross sections for deep inelastic scattering 
and Drell-Yan. The result is then used to write the MS quark coefficient functions in a 
form in which all terms of 0(ln m N) are resummed. 
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Quark initiated bare cross sections, whether for deep inelastic or Drell-Yan processes, 
generally contain logarithmic singularities as z — ► 1, where z is the longitudinal momen- 
tum fraction of the participating partons. These logarithms result from soft and collinear 
gluon emission, and in Mellin space result in logarithms of the form ln m N, which can be 
shown to exponentiate to all orders in perturbation theory Ji]||. Much progress has been 
made [0-0 in resumming these logarithms. Once the general form of the exponent is 
determined, individual coefficients are fixed through matching to fixed order perturbation 
theory. 

Before comparing to large z data it is first necessary to resum the logarithms in the 
factorized quark initiated cross sections, the quark coefficient functions. The resummation 
is then expected to improve the convergence of the perturbation series at large N. In order 
to factorize the cross section it is generally assumed (see for example ref. 0) that in MS 
factorization no resummation is necessary for the quark anomalous dimensions, i.e. that 
higher order contributions to the anomalous dimensions grow no faster than the leading 
order contribution at large N. Considerable support for this assumption was given in ref. 
||, where the large N behaviour of the MS anomalous dimensions was determined to all 
orders using eikonal techniques. 

In this letter we will show that the all order exponentiation ]l],|2j and MS factorization 
are by themselves sufficient to determine the all order large N behaviour of the quark 
anomalous dimensions. More specifically, we take the large N resummed form of the quark 
initiated bare cross sections obtained in ref. 0, then apply dimensional regularization in 
order to factorize the infrared singularities in MS scheme (ref. || uses an explicit infrared 
cutoff). Matching the bare and factorized cross sections then allows us to simply read 
off the asymptotic behaviour of the quark anomalous dimensions: we find that the MS 
nonsinglet qq anomalous dimension ^ q s behaves like O(lniV) at large N, whereas 7^ 
behaves like 0(1/N). Finally, we obtain a general form for the MS quark coefficient 
functions at large N to all orders. The results are then illustrated by explicit computation 
at LL and NLL. 

1. We consider DIS and Drell-Yan bare quark initiated cross sections, which after factor- 
ing off all electroweak factors and decomposing into Lorentz invariants, may be written 
as functions of two kinematic variables, Q 2 and z. For DIS Q 2 is the virtuality of the 
incoming vector boson, while in DY it is the mass of the outgoing dileptons. Similarly, the 
longitudinal momentum fraction z = x = Q 2 /2p.q for DIS, while in DY z = X\X<i = Q 2 /s. 
In both cases we will assume all quarks are massless, and use a s = a s /2ir as a perturbative 
expansion parameter, ignoring all contributions which are suppressed by powers of Q 2 . 



Now soft and collinear gluons radiated from the incoming quark lines can generate terms 
— — '- , 1 < m < 2n, which diverge as z — > 1, and which thus need 
to be resummed if the perturbative cross section is to be improved at large z. In Mellin 
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space, these terms become 
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so the singularities as z — > 1 become logarithmic singularities as N — > oo. Note that in the 
remainder term we do not distinguish terms of O ( ^ n N N ^j from terms of O C^j . In z space 
there are also terms in the bare cross section from virtual graphs of the form 5(1 — z), 
which in iV space are constant at large N. However all the diagrams in which the initial 
quark (or at least one of the two initial quarks in the case of Drell-Yan) is not connected 
to the electroweak boson via a single quark line, i.e. the diagrams contributing to a bare 
pure singlet cross section, are either non singular or of the form ln m (l — z) as z — > 1, and 
are thus of 0(1 /N) after Mellin transformation. Here we will thus be concerned only with 
the nonsinglet diagrams, since it is only these which can have large N singularities of the 
form in eqn. (|l|). Furthermore, we will only be concerned with the singular parts of these 
diagrams (or more precisely those parts which at large N are of 0(ln m N), m = 0, oo), 
which are independent of the type or polarization of the electroweak boson (and thus for 
example in DIS will be identical for all three nonsinglet structure functions -^1,2,3, and in 
DY for dimuon or W production, after removing electroweak factors). 

Consequently, here we need consider only two partonic cross sections (z, Q 2 / k 2 , a s (Q 2 )), 
where the only index a denotes the number of initial state quark lines: a = 1 for DIS, and 
a = 2 for DY. Both partonic cross sections are normalized such that crj^ = 1 for a s = 0. 
Collinear singularities are regulated by a generic infrared cutoff k. In general |],|8[], the 
perturbation series for crj^ in iV space is then (given eqn. ([[])) of the form 

°f (n, %, a s (Q 2 )) = 1 + £ <(Q 2 ) £ C W (Q 2 /* 2 ) In- N + O (1) , (2) 

\ K J n=l m=0 V-iV / 

for some coefficients 6°^ m (Q 2 / k 2 ). Furthermore, the 0(ln m N) terms, for m — 0, 00, not 
only factorize: in Mellin space they also exponentiate, that is the bare cross section may 
be written in the form 

lna^(iV, Q 2 /k 2 , a s (Q 2 )) = (t> [ l\{a s {Q 2 ) IniV, Q 2 /k 2 ) InN 

+ Y,MQ 2 )T<f> [ : ] MQ 2 ) IniV, Q 2 /k 2 ) + O - ) . (3) 

n=0 ViV/ 

This exponentiation provides a convenient framework for organising the expansion: the 
functions (f)^\ contains all the leading logarithms (LL), the next-to- leading logarithms 
(NLL), since there is an extra power of a s , and so on. It is easy to see that in eqn. (||]) 
the LL are then those terms with n + 1 < m < 2n, NLL, those with m = n, etc. Terms of 
O (j^j in will be systematically discarded in what follows, since they are not included 
in its definition. 
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To complete the resummation, we must determine the functions 0) . This may be achieved 
either by eikonal arguments [|1],0|§, or by applying the renormalization group directly to 
the factorization of the large N singularities |J. The result is rather simple: the resum- 
mation may be performed in closed form essentially through a change in the argument of 
the running coupling from a s {Q 2 ) to a s ((l — z) a Q 2 ). Specifically it is found that (see for 
example ref. ||) 



In crM 



V 2V-1 



dz- 



- 1 



1 - z 



J -^A(a s (q 2 )) + B^(a s ((l-z)Q 2 )) 
+K^(a s (Q 2 )), 



(4) 



where the MS renormalization scheme has been used, and soft singularities are regulated 
by the infrared cutoff k. The functions A(a s ), B^(a s ), R^(a s ) each have perturbative 
expansions beginning at 0(a s ). For DY processes the second term is absent, so B^(a s ) = 
0: in ref. || B^(a s ) is simply denoted by B(a s ). The remainder terms K^(a s ) contain 
contributions not necessarily related to soft or collinear gluons. To resum LL and NLL 
logarithms only the first two coefficients in the expansion of A, and the first coefficient 
in the expansion of B^ are necessary, and these may be read off by making a direct 
comparison to the usual fixed order LO and NLO cross sections expanded at large N. 

In order to express the resummation eqn. (|4]) in a form more amenable to subsequent 
discussion, we first rewrite them in d = 4 — 2e dimensions. We may remove the collinear 
regulator by taking k — > 0, since bare cross sections are non singular in a non integer 
number of dimensions. This gives 



by of 
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I %A(a s (q 2 , 6), e) + flW(a,((l - z)Q 2 , e), e] 



+K^(a s (Q 2 ),e) + 0(e), (5) 



for the DIS and DY cross sections respectively. Here a s (fi 2 ,e) is the MS renormalised 
coupling in 4 — 2e dimensions: it satisfies the renormalization group equation 



da s (n 



2 -ea s (/Ae)-£/3 n < + V,4 (6) 

111 ^ 71=0 



where the /3-function coefficients j3 n are independent of e. Clearly a s (/x 2 , 0) = a s (/i 2 ), while 
a s (0, e) = (by analytic continuation from e < 0). The functions A(a s , e), B^(a s , e) and 
K^(a s ,e) all have perturbative expansions in powers of a s , but now the coefficients in 
these expansions will depend on e. The terms of 0(e) in eqn. (|5|) are inconsequential and 
can be dropped. 
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Now using eqn. (|6|) for the running of a s (q 2 ,e), we can do the q 2 integrals term in an 
expansion in powers of a s : 



I 2 a s (Q 2 ,6) 

d\nq 2 (a s (q 2 ,e)) n = J da s 





(7) 



Then from eqn. (|J), we find that the quark initiated bare cross sections at large N are of 
the form 



In a [ ; ] (N,a s {Q 2 , e),e)= / dz 



JV-l 



l-z 



oo 

+ T,9 [ r ] (e)(a s (Q 2 ,e)Y + 0(e). 
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As advertised above, we can now see explicitly the essential feature of the large N resum- 
mation: the change in the argument of the running coupling from Q 2 to (1 — z) a Q 2 . 

2. In this section we will factorize eqn. (|8|) in the MS scheme. The factorization procedure 
involves separating out the collinear singularities from the bare cross sections cxj^ into a 
universal singular factor T q : 



lnal a] = In Cl a] + a In T . 



(9) 



The Cj°] are process dependent coefficient functions, non singular as e — > 0. Just as the 
complete bare partonic cross sections for nonsinglet, singlet and valence processes are all 
proportional to cr' a ' up to terms of 0(1/N), so the complete nonsinglet, singlet and valence 
coefficient functions will all be proportional to up to terms of 0(1/N). Of course the 
same holds true for the singular factor T q , and in particular for the anomalous dimensions 
from which it is constructed: in MS T g satisfies the renormalization group equation 



d\nT q (N } a s Qu 2 ,e),e) 
Sin /i 2 



7 9 (iV, a s (> 2 ,e)), 



(10) 



where 7 g (iV, a s ) = J2^=olq{N)a™ is the anomalous dimension, and in the MS scheme the 
coefficients ^{N) are independent of e. This defines the MS factorization scheme. This 
equation has the usual solution 



T q (N, a s (/i 2 ,e),e) = exp 



d\nq\(N, a s (q 2 ,e)) 



which, when combined with eqn. (0) and expanded in powers of a s generates all the 
collinear singularities of the bare cross section, in the form of inverse powers of e. 
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We first use eqn. (||) to determine the degree of divergence of the coefficients (e) and 
g\ a \e) as e — > 0. Differentiating eqn. (Q) with respect to lnQ 2 , and using eqn. ([X0|) , we 
find 



<91naW 91nCW / q 



(12) 



Since both terms on the right hand side are nonsingular as e — > 0, it follows that 
9 In crj^ /<9 In Q 2 is non singular. Differentiation of eqn. (|8|) then leads to the conclusion 
that 



if(e) =£/S<r* + 0(e), <^(e) =£sge-* + 0(e 



t=o 



(13) 



as well as various relations among the f- a t , and . 
Using eqn. (0), we can rewrite lnl^ in the form 



1 Q rl 2 

\nT q (N,a s (Q 2 ,e),e) = J ' dzz^ 1 J -^P q (z, a s (q\ e)) 



(14) 



where P q (z,a s ) = Y^=iPq{ z ) a% s is the quark splitting function (so the P q (z) Mellin 
transform to the anomalous dimensions 7™). Since the 7™ (N) are independent of e, so too 



are the P q (z). Substituting eqns. (fb§,BI) into eqn. (0) we then find that 



\ r 1 00 i 
\nCf= dz{z N - l -l) ££4?e-<M(l-*rQ 2 ,e)) 1 

-a / -|-P g (z,a s (g 2 ,e)) 
^ 
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/ dz{z N - 1 - 1) — E^ ] (a s ((l - z) a Q 2 ,e)r + E^(a,(Q 2 , e)) 

-a j - 1) j %P q {z, a s (q 2 , e)) 



+ 



+ 



(l-z) a Q 2 



+ °(^J+°( e )- (15) 



The motivation for organisation of the terms in this equation is that, as we will now show, 
each of the three pairs of square brackets is finite as e — > 0. 



5 



We first work with the first pair of square brackets. Using eqn. (|6]), we can perform the 
integration of the a™(q 2 ) in the second term as 



(l-z)"Q 2 2 a s ((l- 2 )°QV) 

— (a s (q 2 ,e)) n = / da s ^ 

oo 



£ Kn (e)(a s ((l-zrQ 2 ,e))\ (16) 



It follows that 
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We now expand eqn. ( fP7D in e. In eqn. (0), the Pi, n (e) may be expanded as 

i—n+l 

PiM = E Ptnt~ S - (18) 
s=l 

Note that p s in — if s = 1 and i ^ n, but we do not explicitly show this in order to 
simplify the notation. Substitution in eqn. ( |T7D then gives finally 

(1— z) a Q 2 ^ oo i i—n+l 

[ -^P^,a s (y 2 ,e))=$:E E W^-> s ((l-*) QV)r 

q y i=l n.=l s=l 

oo i i— i+1 

= E E E TOp^ - ' (a.((l - ^) a Q 2 , e))\ (19) 

i=l t=l n=l 

The second term in the second pair of square brackets in eqn. ([15]) may be expanded 
similarly, using eqn. (|19| ) with a = 0: 



Q dg 2 x ' 



/ ^P g (z, a s (g 2 , e))=EE E W^MQ 2 , ^)) 1 - (20) 

n y i=l t=l n=l 



Finally, the last term in the third pair of brackets in eqn. (15) is given by the difference 
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of eqns. ( |I9|) and (p0|), and is thus finite as e — > 0. Indeed, we can write it as 

q 2 <3 2 2 oo 

/ ^P g (^a s (g 2 ,e))= / ^ EK^))" + O(e). 



(!-*)■ 



(l-^) a Q 2 



Putting all this together, substituting eqns. (|i~9| pOlpll) into eqn. (p~5f ) we find 
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+ 0(e). 



(22) 



Now, consider only the coefficients of e~* for t = 1, oo in eqn. ([22|), which must vanish 
so that is finite in the limit e — > 0. In the first pair of square brackets, these coefficients 
are all functions of N, in the second pair of square brackets they are all independent of N, 
while all terms in the third pair of square brackets are nonsingular. Thus each of the three 
pairs of square brackets is separately nonsingular as e — > 0. Moreover the cancellation of 
singularities in the first pair of square brackets implies that, for i > t > 1 and z ^ 1, 



1 Z n=l 



(23) 



We can treat p\ n in eqn. ( |23"D as a set of matrices with indices % and n, but vanishing for 
n > i — t + 1; each matrix is then triangular and thus invertible. Then multiplication of 
both sides of eqn. (|23|) by the inverse of this matrix gives an expression with just P™ on 
the right hand side, with the left hand side proportional to (1/(1 — z)). Thus we find that 
for 2 / 1 



1 



P q (z,a s ) = Q(a s ) — 



(24) 



where Q{a s ) = YhZi a l s Qi, and the Qi may be determined by substitution into eqn. fl23]) . 
Note that for consistency we also require that f\ a } is proportional to a for t > 0: this is the 
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essence of universal factorization. Similarly the cancellation of singularities in the second 
bracket implies that 



i-t+l 1 



$ = a E / dzP?(z)pl n , (25) 

n=l i 



and thus that Jq dzP™{z) is finite. It follows that for all z, we may write 



1 

1 



P q (z,a s ) = Q(a s ) 



l-z 



+ R{a s )S(l-z), (26) 



where again R(a s ) = a l s Ri, and the Ri may be determined by substituting eqn. (|26|) 
into eqn. (|25|) and inverting. 

Thus we find that, in Mellin space, as iV — * oo, the anomalous dimension 

lq {N, a,) = -Q(a s )(lnN + lE ) + R(a s ) + 0(1/N), (27) 

where je is Euler's constant. Remembering that all the nonsinglet, singlet and valence 
anomalous dimensions are equal to 7 9 up to terms which vanish as 1/N in the large N 
limit, then since in the usual decomposition 

=S-^ NS + ^ PS ry _ =$..ryNS + PS NS = NS ± NS (2g) 

I lilj U V Iqq ' Iqq ' IQiQj Ul J Iqq ' Iqq ' '± Iqq > qq ' V °/ 

we must have 



lqq{N) = lq (N) + 0(1/N), ^(N) = 0(1/N), 

fqq S {N)=0{l/N), 7 « 



< 8 {N) = 0(1/N), Tff{N) = 0(1/N), (29) 



consistent with our remarks earlier that diagrams in which the quark evolves into a gluon 
are suppressed by 1/N. 

It remains to take the limit e — > in eqn. (|22]) to give an explicit expression for the 
resummed large iV MS coefficient function: we find 



In Cl a 



JV-l 



dz- 



1 r 



{l-z)" 



l-Z 



dq 



2 oo 



E Qn(a s (q 2 )) n + E fl>s((l ~ z) a Q 2 )Y 



n=l 



i=l 



N 



(30) 
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We immediately recognise this as eqn. (£|), with 



A(a s 



oo 

E 

n=l 



Qi.a 



B [a] 



E4o< 



oo 

i=i 



(31) 



and with the infrared regulator set equal to the factorization scale, i.e. k 2 = Q 2 . 

To calculate the coefficients /Jo and the gj^j, we simply expand eqn. ( P0| ) in a s (Q 2 ), i.e. 
undo the resummation, and compare coefficients of a l s (Q 2 ) with those in the fixed order 
coefficient functions. The universal coefficients Qi of the O(lniV) terms in the anomalous 
dimension serve as a consistency check. This procedure will be illustrated explicitly for 
all LL and NLL terms in the next section. 

3. We now show explicitly that the factorized coefficient functions eqn. (|3Qf) correctly 
resum all LL and NLL terms, by comparing then to the fixed order NLO and NNLO 
coefficient functions and thereby deduce the leading behaviour at large N of the LO and 
NLO quark anomalous dimensions. 

To NLO, eqn. © reads 



In Cl a \N,a s (Q 2 )) 



JV-l 



dz- 



aQ l ^fa s ((l-z) a Q 2 ) 



a s (Q' 



+a - (a.((l - zYQ 2 ) - a s (Q 2 )) + /ga.((l - z) a Q 2 ) 

+0(a s (a s \nN) m ). 



(32) 



Performing the z integral, we find 



\nCW(N,a s 




ln(l — a\ s ) + a\ s 

A) 

- a\ s ) 




- a\ s ) 

aX s + 0{a s {a s In iV) m ), (33) 



where X s = a s /3 lniV. This result should be compared to the general expression eqn. (g): 
the first line gives the LL terms 0_x, while the rest gives the NLL </> . 

To determine the large N behaviour of the LO and NLO anomalous dimensions, i.e. the 
coefficients Q\ and Q2, and also the nonuniversal coefficients /j"d, it is sufficient to expand 
eqn. (^) in a s (Q 2 ) to NNLO and compare to NLO and NNLO coefficient functions at 
large N. The result, after exponentiating, is 
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Cf (N,a,) = l + a s Qa 2 Q x In 2 N + (a 2 Q llE - /fj) lniv) 
+a 2 (^-a 4 Q 2 In 4 iV + [^a 3 g x /3 + +a 2 Qi (a 2 Q l7£ - /$ 



In 3 AT 



+ 



In 2 N) + 0(a 3 ). (34) 



Now, for large JV, to NLO [ 10 1 and NNLO Oj the DIS quark coefficient function 



Cl 1] (N,a s 



1 + a 

y2 



C F In 2 N + C F (2 lE + |) IniV + 7 2 + ^ 7e - C(2) + O (1 
3„ 2 11 



2C^ 7i; - -C^n/ + ~C£ + — C^C^J ln d iV 



+ (yCVC A7 £ + ^ 2 7£ - C7 J? (7 A C(2) - yC| - C 2 C(2) + ^CVC A 

2 29 \ 

+3C| 7 | - -C F T R n flE - —C F T R n f ) ln 2 iV + O(lniV) 



(35) 



Comparing this with eqn. ( p4l) with a = 1, we find 

'67 



CJi — 2C F , Q 2 — C f Ca 



2C(2) 



20 



C F T R n f , 



(36) 



consistent with the large iV behaviour of the LO [E| and NLO |13| anomalous dimensions. 
Furthermore the coefficient 



Ji,o — — 2 



(37) 



as in @J|. For Drell-Yan at large JV, to NLO |J and NNLO |5 



Cf{N,a, 



1 + a s 



4C F In 2 iV + 8C f7e In N + ^C Fl 2 E + 8C F C(2) - 8C F + 0(1) 



/44 16 
8C 2 In 4 JV + [-^C F C A - —C F T R n f + 32C 2 7b ) In 3 N 



40 16 44 134 

- —C F T R n f - —C F T R n flE + —C F C AlE + — C F C A - 4C F C A ((2) 

+32Cf<(2) + 48C| 7 | - 32C 2 F ) In 2 N + 0(lnN) 



(3? 



which is again consistent with eqn. (|34|) , this time with a = 2, provided eqns. ( p6|) for the 
large A" anomalous dimensions are satisfied, and /] 2 q = 0. With this result and eqn. fl3~7|), 
eqn. (|33"D agrees with the results of 0] at NLL. Moreover, our general proof of eqn. (|26|) 
now places the NNLL results of on a firmer footing. 

4. We have used the dimensionally regularized form of the large A" resummed bare quark 
initiated cross sections for DIS and Drell-Yan, eqn. together with the fact that the 
bare quark initiated cross sections are independent of the species of the electroweak boson, 
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and the definition of MS factorization, to obtain to all orders the O(lniV) behaviour of 
the nonsinglet MS anomalous dimensions P^ q s and the O (jf) behaviour of P^ s and the 
pure singlet. This has interesting implications for the large x evolution of MS quark parton 
distribution functions: at large N the evolution factor 

r(Q 2 //i 2 ) = N'^ 2 ~k Q( - a ^\l + 0(l/N)). (39) 

so that if q(x,fi 2 ) ~ (1 — x) b ^ 2 ^ as x — > 1, this behaviour persists at higher scales with 

Q 2 , 2 

b{Q 2 ) = % 2 ) + f -\Q{a s {q 2 )) (40) 
p? 

order by order in perturbation theory. 

We then used the result (^6]) to construct a general large N resummed expression for the 
MS quark coefficient functions for DIS and DY, eqn. (p0[), and in particular showed that 
all the large N singularities in the DY case can be deduced from the O(lniV) terms in 
the quark anomalous dimension. We verified these results explicitly at LL and NLL. For 
the future, the large N behaviour of NNLO anomalous dimensions and NNNLO DIS and 
DY coefficient functions will provide a useful consistency check on new calculations. 

Acknowledgements: SDA would like to thank G. Sterman for very useful discussions 
on this topic, and W. van Neerven and S. Catani for correspondence. 
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